final_full <- final %>% mutate(days_scaled = as.numeric(scale(days, center = F)))
m1 <- lm(violence ~ days_scaled * treatment, final_full, weights = weight)

m1a <- lm(violence ~ days_scaled * treatment+age+race6+educ4+gender, final_full, weights = weight)

m2 <- lm(norms ~ days_scaled * treatment, final_full, weights = weight)

m2a <- lm(norms ~ days_scaled * treatment+age+race6+educ4+gender, final_full, weights = weight)

m3 <- lm(affpol ~ days_scaled * treatment, final_full, weights = weight)

m3a <- lm(affpol ~ days_scaled * treatment+age+race6+educ4+gender, final_full, weights = weight)

stargazer::stargazer(m1,m1a,m2,m2a,m3,m3a, type = 'latex',
                     title="Full Models",
                     omit.stat = c("ser", "adj.rsq", 'rsq'), 
                     font.size = "scriptsize", column.sep.width = "6pt",
                     df = F, no.space = T,
                     dep.var.labels = c(
                       "Violence", 
                       "Norms", 
                       "Aff Pol"))
